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Abstract 

We construct numerically static non-uniform black string solutions in 
six dimensions by using pseudo-spectral methods. An appropriately de¬ 
signed adaptation of the methods in regard of the specific behaviour of 
the field quantities in the vicinity of our numerical boundaries provides 
us with extremely accurate results, that allows us to get solutions with an 
unprecedented deformation of the black string horizon. Consequently, we 
are able to investigate in detail a critical regime within a suitable param¬ 
eter diagram. In particular, we observe a clearly pronounced maximum in 
the mass curve, which is in accordance with the results of Kleihaus, Kunz 
and Radu from 2006. Interestingly, by looking at extremely distorted 
black strings, we find two further turning points of the mass, resulting in 
a spiral curve in the black string’s phase diagram. 
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1 Introduction and Summary 

When Gregory and Laflamme (GL) discovered the instability of uniform black 
strings m (UBS), a search began for a new branch of solutions that emanates 
from this instability. Such non-uniform black strings (NBS) were constructed 
first perturbatively [aisiii] and later numerically 0151012] by several authors. 
In Ref. [2] a natural measure of the deformation of the black string horizon, 
^ = 5 '^mhi('^max — ^min), was introduced, where i?max (.Rmin) is the maximum 
(minimum) of the NBS’s areal horizon radius along the compact dimension (for 
A = 0 we thus have the UBS). 

An interesting issue in this context is the conjectured phase transition of 
the NBS leading to a localized black hole [8] (BH), occurring when the black 
string’s horizon has a critical deformation (A —^ oo) and pinches off. So far the 
numerics gave strong evidence in favour of this conjecture [nmniiii], although 
the actual transition point was not reached, see Refs. [UlSj. 

The aim of this work is the construction of highly accurate numerical solu¬ 
tions in order to clarify the controversy about the existence of a maximum in the 
NBS mass curve (cf. Refs. 00 [ 7 ]), at least in six dimensions. For this purpose 
a pseudo-spectral scheme was implemented, with sophisticated adaptations to 
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get satisfactory accuracies in the large A regime. The techniques include the use 
of several appropriate coordinate mappings, the introduction of multiple do¬ 
mains and the split of each metric function into two parts (near infinity). Also, 
a large number of grid points near the critical point (located on the horizon) 
was to be taken. Consequently, we were able to reach values of A ~ 200. 

The described method provided us with results that not only support those of 
Ref. [3 and therefore confirm the maximum in the mass curve. Moreover, they 
show two further turning points in this curve, i.e. in sum we found three turning 
points for the mass, and likewise for the relative tension, the temperature and 
the entropy. Consequently, we obtain the beginning of a spiral curve in the NBS 
phase diagram for large A. Since our method ceases to converge for values of 
A > 200, possible further turning points are to be discussed elsewhere. 

The present results are restricted to six dimensions. Note however that 
work is in progress regarding five dimensions. The corresponding problem is 
numerically much more difficult to deal with since the metric functions show 
logarithmic type behaviour, which becomes a subtle issue if high accuracy is 
needed. Another interesting question, to be treated in a forthcoming article, 
regards the occurrence of a spiral curve in the BH phase diagram, which shows 
a continuous transition to the NBS spiral. If this occurrence can be confirmed, 
then the phase transition conjecture would have further evidence. Furthermore, 
we note that a search for unstable modes corresponding to the turning points 
in the phase diagram would be illuminating. 


2 Metric and Charges 


We consider the static NBS metric in D dimensions and with the background 
^ gi form 


ds2 = _e2^(’''^)/(r)dt2 (1) 

The three unknown metric functions A, B and C depend on the radial coordinate 
r S [ro,oo] and the coordinate z S [0,T] varying along With /(r) = 1 — 
(ro/r)^“^, the horizon resides at r = rg. All functions are periodic in z with 
period L. Additionally, the functions A, B and C possess reflection symmetry 
in z with respect to the coordinate value z = L/2. 

For A = B = C = d, Eq. ([1]) describes the UBS. Then, the GL instability 
occurs, if the size of the compact dimension L exceeds a certain value Tgl- 
The mass M and relative tension n are obtained from the leading coeffi¬ 
cients of the asymptotic behaviour of A and B when r oo. In contrast, the 
temperature T and the entropy S can be read off from horizon values (see, for 
instance. Ref. 0)- These four quantities obey Smarr’s formula |14) 


as well as the first law of thermodynamics 

dM = TdS + ^dL. 

Lj 


( 2 ) 

( 3 ) 


In the following we consider D = 6. Here the critical value of the size of the 
compact dimension is Lgl = 4.95161542007(1) • rg (see also footnoted]). 
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3 Numerical Method 


The basis of our numerical scheme is the expansion of each function in terms 
of Chebyshev polynomials. Considering the functions’ values on Lobatto grid 
points (which include the boundaries), we can easily calculate spectral deriva¬ 
tives. To solve a system of differential equations, a Newton-Raphson method is 
applied. 

From Einstein’s vacuum field equations one finds the system of partial dif¬ 
ferential equations together with some boundary conditions (see Ref. [J] for a 
detailed discussion). A unique solution can be obtained by scaling physical 
quantities in terms of appropriate powers of tq and fixing L (we put L = Tgl), 
and by prescribing a value of the function B on the horizon at z = Ll2. 

By an analysis of linear perturbations around the UB^I we find that the 
following ansatz leads to a rapid decay of the spectral coefficients: 
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According to this ansatz, we have to solve for three functions Ai, Bi and Ci, 
which depend on the two coordinates r and z. Moreover, the functions Aq, 
Bq and Co, depending merely on the radial coordinate r, are to be found. A 
particular benefit of the decomposition Q is the fact that M and n follow from 
the asymptotics of Aq and Bq alone. 

In addition to the aforedescribed split, we introduce new coordinates ^ and 
u via the transformations rg/r = ^^(2 — ^)^ and u = cos(27rz/L), in which the 
horizon is located at ^ = 1 and infinity at ^ = 0, and where z = Lj^ corresponds 
to M = —1 and z = 0 to M = 1. Note that the term (rg/r)^/^ appearing in ® 
is regular with respect to Also, all ^-derivatives of the occurring functions 
vanish, when taken at the horizon ^ = 1. 

With this ansat^ we are able to produce accurate solutions in the NBS 
branch. However, as in the limit of infinitely pronounced horizon deformations 
(i.e. A oo) the values of the functions diverge at the critical point (^, u) = 
(1, —1), we introduce the auxiliary functions a = exp(—2A), /3 = exp(—2i?) and 
7 = exp(2C'). They stay finite within the entire integration domain. Neverthe¬ 
less, steep gradients appear in vicinity of the critical point, for which reason we 
introduce appropriate coordinate mappings in order to increase the resolution 
near this point[l, see Fig. [T] 

^We set A = ea, B = eb, C = ec for some small e and considered the first order in e. 
The analysis of the linear problem also reveals the value of Lql (see appendix in Ref. US]) 
and provides an initial guess for the Newton-Raphson scheme in the non-linear regime. We 
determined Lql units of ro) up to a precision of 10“^^, cf. section[2l 

^Note that we utilized the additional split Co(0 = ^00 +^^C'oi(^) which proved to be 
useful. 

^We note that the coordinate u is not suitable in this region for large A (u-derivatives 
become very large at the critical point). We cured this issue by going back to the coordinate 
2 :, or, to be more precise, to a rescaled version of \t, u = 1 — ^z. 


3 




Figure 1: We divide the integration domain up into seven sub-domains. For ^ < 0.5 
we solve the equations with the ansatz Here, the three sub-domains guarantee a 
fast convergence of the spectral coefficients in each of these domains (this subdivision 
helped to deal with the non-analytic exponential terms appearing in which occur 
in higher orders in the non-linear regime). Another benefit is that the resolution 
in M-direction can be adapted in each domain (remember that the u-dependency is 
exponentially suppressed for ^ ^ 0 , while for ^ far from zero high resolution in u 
is needed). For ^ > 0.5 we solve for a = e~^^, fi = e“^® and 7 = e^*^ (see text 
for explanation) and introduce further coordinate mappings, that lead to convergence 
of the coordinate lines in the vicinity of the critical point. An enhancement of the 
resolution is achieved by a further division into sub-domains. 


4 Results 

With the methods described above we are able to construct solutions with the 
deformation of the black string horizon reaching up to A « 200. Note that for 
large A, the accuracy obtained is still of the order 10“^^ (absolute) and 10“® 
(relative). As an additional accuracy check we evaluate Smarr’s relation and 
found agreement of left and right hand sides of ([5]) to similar orders (see Fig.®). 
The first law dM = TdS (remember dL = 0, since L was fixed) is also satisfied 
to great precision over the whole range of solutions. 


Am 200 



N 

Figure 2 : Convergence of the residuum Rn and the deviation from Smarr’s formula 
Asmarr for the solutiou with the highest deformation achieved, A « 200. We calculate 
Rn hy comparing a solution with high resolution to a solution with smaller resolution 
N on a fine grid and take the maximal difference. 
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The high accuracy achieved permits the detailed consideration of particular 
black strings’ phase diagrams, see Fig. [31 Our results are in very good agreement 
with those of Ref. , which show a maximum in the mass curve (as well as in 
the entropy curve, whereas relative tension and temperature respectively show 
a minimum). Going beyond these extremal points by increasing the degree 
of deformation, another two turning points appear in each of these diagrams. 
We thus encounter the beginning of a spiral curve in the black strings phase 
diagrams. Particularly, in Fig. [3] a spiral with about one and a half turns in the 
(n, S) and (T, M) diagram can be seen (though it becomes tiny). It is tempting 
to speculate that there are infinitely many turns when going to A ^ oo, which 
would indicate a growing number of unstable modes. 
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Figure 3: Left: Entropy as a function of relative tension. Right: Mass as a function 
of temperature. The top panels show all of our data. From top to bottom, the panels 
show more and more magnified portions of the spiral in question. All values are 
normalized with respect to their corresponding values of the UBS (denoted with lower 
index 0). 
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